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ABSTRACT 


An  inviscid  theory  is  developed  for  the  steady  aerodynamic  loading  on  a 
ducted  propeller  in  forward  flight  at  angle  of  attack.  As  a  first  approximation, 
the  problem  is  regarded  as  the  superposition  of  the  given  ducted  propeller  at 
zero  incidence  and  a  cylindrical  duct  of  the  same  aspect  ratio  but  of  zero 
thickness  and  camber  at  the  given  incidence.  This  approximation  is  then  im¬ 
proved  to  include  the  cyclic  variation  of  the  blade  loading.  Net  only  are  the 
normal  force  and  pitching  moment  changed  but  also  a  finite  side  force  and  yawing 
moment,  unknown  heretofore,  are  found.  Explicit  formulas  are  given  to  carry  out 
numerical  calculations  and  a  limited  comparison  with  experimental  results  show 
good  agreement. 
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PRINCIPAL  NOMENCLATURE 


unsteady  airfoil  coefficients;  see  Eqs.  (45)  and  (58)  and  Table  11 


a  ,b 
v'v 


Glauert  coefficients  for  sin  9  and  cos  9  portions,  respectively, 
of  steady  shroud  vortex  loading;  see  Eqs.  (4),  (49)  and  (50) 


sin  9  and  cos  9  portions  of  steady  shroud  vortex  loading  respec¬ 
tively;  see  Eq.  (4) 


shroud  chord 


shroud  surface  pressure  coefficient;  see  Eq.  (24) 


axisymmetric  portion  of  cp  ,  identical  to  cp  for  zero  incidence 

Glauert  coefficients  for  axisymmetric  portion  of  steady  shroud 
vortex  loading  7q/U 

shroud  and  propeller  force  and  moment  coefficients;  see  Eqs.  (31) 
and  (32),  for  example 


blade  chord 


£  shroud  source  strength  per  unit  axial  length 

«  Fourier-type  coefficients  of  {jf/ d  and  ®/a  respectively;  nee 
rv  v  Eqs.  (51)  and  (52) 

Si  sin  9  and  cos  9  portions,  respectively,  of  steady  radial  velocity 

induced  on  shroud  reference  cylinder  due  to  cyclic  loading  on  pro¬ 
peller;  see  Eqs.  (46)  and  (48; 

J  propeller  advance  ratio,  u/fiR 


reduced  frequency  of  sinusoidal  inflow  to  propeller  blade  sections; 
see  Eq.  (44) 


kernel  of  integral  equations  for  inclined  duct  problem;  see  Eqs*  (6) 
and  (19) 

blade  index;  l  ■  0,  1,  ...  ,  (N-  1) 
shroud  lift 


propeller  lift 
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u,v,w 

AAA 

u,v,w 
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vf  ,7#  ... 
Vf0'70 


shroud  pitching  moment,  positive  nose-up 

propeller  pitching  moment,  positive  nose-up 

shroud  yawing  moment,  positive  nose-right 

number  of  propeller  blades 

propeller  yawing  moment,  positive  nose-right 

inversion  matrix;  see  Bqs.  (23),  (53)  and  (5*0,  and  Tables  2-5 

static  pressure  on  shroud 

static  pressure  at  infinity 

Legendre  function  of  second  kind  and  half  order  (n+*s)  with 
argument  a 

reference  radius,  taken  as  mean  shroud  radius  at  propeller  plane 
propeller  radius 

shroud  aide  force,  positive  in  positive  y-direction 
propeller  side  force,  positive  in  positive  y-direction 
extended  Hough  matrix;  see  Eq.  (27)  and  Tables  6-9 
time 

shroud  thickness  distribution 

steady  x,r, 6 -components  of  disturbance  flow  respectively 
unsteady  x,r, 6 -components  of  disturbance  flow  respectively 
flight  speed 

steady  portions  of  radial  velocity  induced  by  f,7,  ...  respectively 

axisymmetric  portions  of  v£  ^  respectively 

resultant  velocity  (02+n2r2),J  relative  to  blade  element 
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x,x,e 


x,y,z 


*P 


x,r 
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V7X 
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V*w 


a'°b 
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(<p}«  (0) 
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shroud-fixed  cylindrical  coordinates;  see  Pig.  1 
shroud-fixed  Cartesian  coordinates;  see  Pig.  1 

axial  propeller  position,  positive  downstream  of  shroud  midchord; 
see  Fig.  1 

x/R  and  r/R  respectively 

angle  of  attack  to  free  stream,  in  radians 

total  steady  bound  shroud  vortex  strength  per  unit  axial  length 

axisymrcetric  and  S -dependent  portions  of  y  respectively;  see  Eg.  (4) 

steady  and  unsteady  portions,  respectively,  of  bound  circulation  dis¬ 
tribution  on  fth  wlade 

local  inclination  angle  of  shroud  camber  line  relative  to  x-axls; 
see  Pig.  1 

angular  por  tion  of  1th  blade,  (-ftt+  Shri/N) 

phase  angles  for  unsteady  blade  circulation  and  lift  per  unit  length 
respectively;  see  Eqs.  (41)  and  (43)  and  p.  27 

shroud  chord-to-diameter  ratio,  c/2R 

V" 

inflow  velocity  factors;  see  Eqs.  (36)-  (38)  and  Table  10 

amplitude  factors  for  unsteady  blade  circulation  and  lift  per  unit 
length  respectively;  see  Eqs.  (41)  and  (42)  and  p.  27 

chordwise  Glauert  variable;  x  ■  -  Xcos  0 

column  vectors  associated  with  three-dimensional  contribution  to 
Cp  ;  see  Eqs.  (27)-  (29)  and  (55) 

angular  rotational  speed  of  propeller,  radians  per  unit  time 
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INTRODUCTION 


In  previous  studies  wa  have  treated  the  ducted  propeller  extensively  for 

l_ii 

the  case  of  forward  flight  at  zero  incidence  ,  as  well  as  static  and  low- 
speed  flight^.  Sere  our  attention  is  turned  to  the  case  of  forward  flight  at 
angle  of  attack,  wo  take  the  design  propeller  loading  and  the  shroud  geometry 
as  known  and  seek  to  determine  the  steady  part  of  the  shroud  loading. 

Consider  first  the  related  problem  of  the  isolated  duct  at  angle  of  attack, 

/f  7  A 

which  has  been  studied  in  detail  by  B.  Ribner  ,  J.  Weissinger' '  and  many  others. 
From  Weissinger 1 s  theory,  the  duct  at  incidence  may  be  regarded  as  the  super¬ 
position  of  the  given  duct  at  zero  incidence,  plus  a  cylinder  of  the  same  aspect 

ratio  but  of  zero  thickness  and  camber,  at  the  given  incidence.  The  results 

q  in 

have  been  found  to  be  in  excellent  agreement  with  experiment*7' 

The  ducted  propeller  at  angle  of  attack,  however,  is  considerably  more 

difficult.  This  is  clear  at  the  outset  on  the  basis  of  the  difficulty  encoun- 

11  12  1^ 

tered  in  dealing  with  the  bare  propeller  at  incidence  '  .  A.  R.  Kriebel  has 

treated  the  problem  by  an  a  priori  extension  of  the  above  superposition  princi¬ 
ple.  That  is,  he  regards  the  ducted  propeller  at  Incidence  as  the  superposition 
of  the  given  ducted  propeller  at  zero  incidence,  plus  a  cylinder  of  the  same 
aspect  ratio  but  of  zero  thickness  and  camber  at  the  given  incidence.  Bis  ap¬ 
proach  does  not  seem  too  unreasonable  and  reduces  the  problem  to  the  sum  of  two 

2i  Q 

problems  for  which  the  solutions  are  available  .  unfortunately,  though,  it 
completely  ignores  the  incidence  of  the  flow  to  the  propeller  plane  and  hence 
the  cyclic  fluctuation  of  the  circulation  on  the  propeller.  For  a  bare  propeller 
this  cyclic  fluctuation  causes  additional  steady  loads.  As  a  result,  we  expect 
similar  loads  for  the  ducted  propeller,  not  only  on  the  propeller  ibielf  but 
also  on  the  duct  due  to  its  interaction  with  the  propeller. 

We  start  with  a  general  formulation  of  the  problem  and  then  examine  the 
superposition  model  as  a  first  approximation.  The  cyclic  fluctuation  of  the 
loading  on  the  propeller  and  the  associated  interaction  with  the  shroud  follow 
and  an  improved  solution  together  with  illustrative  results  are  found.  In  con¬ 
clusion,  comparison  with  seme  DTKB  experiments  are  made. 
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CHAPTER  ONE 


GENERAL  FORMULATION 


1.1  Physical  Model  and  Problem  statement 


We  consider  a  ducted  propeller  in  steady  forward  flight  U  with  a  constant 

rotational  propeller  speed  Cl  and  angle  of  attack  a  in  an  unbounded,  inviscid, 

incompressible  fluid  otherwise  at  rest,  see  Fig.  1  .  The  coordinate  system 

1  2 

(x,r,0)  shown  is  fixed  to  the  shroud.  In  our  previous  work  '  for  the  case  of 
zero  incidence,  propeller-fixed  coordinates  were  used  since  the  flow  was  steady 
with  respect  to  such  a  frame.  In  the  present  case,  however,  the  flow  is  unsteady 
with  respect  to  any  reference  frame  and  shroud-fixed  coordinates  are  found  to  be 
more  convenient. 

We  assume  that  both  the  angle  of  attack  and  the  shroud  camber  are  small. 

We  also  assume  that  the  shroud  thickness  is  small  compared  to  the  shroud  Chord 
c  .  The  a  propeller  blades  are  taken  to  be  slender  and  lightly  lorded,  and 
the  hub  is  neglected.  As  a  reference  radius,  we  choose  the  radius  R  of  the 
camber  line  at  the  propeller  plane. 

With  these  assumptions,  we  may  represent  the  propeller  by  the  classical 
lifting-line  model  Wherein  each  blade  is  replaced  with  a  single  bound  radial 
vortex  and  an  accompanying  helical  sheet  of  trailing  vortices.  The  situation 
is  complicated,  though,  by  the  fact  that  the  blade  circulation  must  experience 
a  cyclic  fluctuation  about  its  design,  or  mean,  value  due  to  the  incidence  of 
the  ducted  propeller  to  the  free  stream.  Consequently,  additional  shed  vor¬ 
tices  enter  the  stream  at,  and  parallel  to,  the  lifting  lines  and  convect 
downstream. 

As  in  thin  annular  airfoil  theory,  a  superposition  of  ring  sources  and 
bound  ring  vortices,  together  with  an  appropriate  system  of  trailing  vortices, 
is  used  for  the  shroud.  Again,  the  situation  is  complicated  by  tho  addition  of 
shed  vortices  since  the  bound  vortex  strength  varies  with  time.  This  arises 
partly  from  the  cyclic  fluctuation  in  the  propeller  loading  and  partly  from  the 
change  in  the  radial  wash  induced  by  the  mean  loading  itself  as  the  blades 
rotate. 

With  the  above  model  in  mind,  the  problem  may  be  stated  mathematically  as 
follows:  Given  the  geometry  of  the  configuration  and  the  mean  impeller  blade 
circulation,  determine  the  shrcud  source  distribution  and  the  steady  part  of  the 
bound  shroud  vortex  distribution  such  that  the  tangent-flow  condition  is  satis¬ 
fied  at  the  reference  cylinder  and  the  Kutta  condition,  at  the  trailing  edge. 

1.2  Tangont-Plow  Condition 

The  linearized  tangent-flow  condition  is  well  known  and  has  the  form  given 
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e  ±  Jj  -  a  sin  0  =  ^  -c/2  <  x$  c/2  ,  r  -  R+  (1) 

where  e  denotes  the  local  inclination  of  the  cumber  line  to  the  x-axis,  and 
t  denotes  the  local  shroud  thickness,  the  prime  indicating  differentiation 
with  respect  to  its  argument  x  .  In  addition,  we  take  u,v,w  and  u,v,w  to 
be  the  steady  and  unsteady  axial,  radial  and  circumferential  portions  of  the 
disturbance  flow  respectively.  Eq.  (1)  is,  in  fact,  only  the  steady  part  of  the 
tangent-flow  condition.  The  unsteady  part  is  simply  0  *  v  on  the  mean  duct 
surface,  similarly,  the  steady  and  unsteady  portions  of  the  flow  must  satisfy 
the  Kutta  condition  separately,  since  the  total  flow  must  satisfy  the  Kutta  con¬ 
dition.  In  seeking  the  steady  shroud  loading  we  shall  be  concerned  only  with 
Eq.  (1)  together  with  its  associated  steady  Kutta  condition. 

1.3  Evaluation  of  the  steady  Radial  Hash 

The  distribution  of  the  shroud  source  and  vortex  strengths  per  unit  length 
sure  taken  to  be  f  and  (7+7)  respectively,  and  the  bound  circulation  dis¬ 
tribution  over  the  ith  blade  is  defined  as  (r+  r, )  .  Of  these,  only  V  is 
specified  a  priori. 

The  trailing  vortices  from  the  shroud  -d (7  +  7 )/dP  are  assumed  to  convect 
downstream  with  the  speed  U  .  Each  bound  ring  vortex  therefore  contributes  a 
semi -infinite  cylinder  of  straight  trailing  vortices.  The  corresponding  shed 
vortices  -dy/dt  constitute  a  similar  cylinder  but  composed  of  ring  vortices. 

The  trailing  vortices  from  the  propeller  -d(r+f^)/dr  are  also  assumed  to 
convect  downstream  with  the  speed  U  and  so  form  regular  helices  with  advance 
ratio  O/flr  .  The  shed  vortices  -df^/dt  emitted  from  each  blade  fall  on  these 
helices  and,  as  above.  Me  perpendicular  locally  to  the  trailing  vortices. 

The  steady  radial  wash  v  may  be  expressed  accordingly  as  the  superposi¬ 
tion  of  the  individual  contributions,  or 

v“  vf+  v7+  v^+  vy' +  v-$'+  v-f,  +  vr+  +  vr'+  vf '+  vf"  (2) 


Where  the  prime  refers  to  the  trailing  vortices  and  the  dot  to  the  shed  vortices; 
vf  is  the  steady  radial  wash  induced  by  the  eource  distribution,  and  so  forth. 

Now  consider  each  contribution  in  turn,  we  know  vf  is  independent  of  9 
and  hence  denote  this  fact  by  setting 


(3) 
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This  quantity  ia  calculated  explicitly  in  Ref.  3  ,  but  for  our  purposes  the 
above  will  ouffice. 

To  evaluate  the  shroud  vortex  contributions,  wo  take  the  steady  part  7 
in  the  form, 


y  ■  y0  *  7J. 

*>  VC  +  UA  sin  0  -fr  HB  cos  0  (4) 


with  c  ,  A  and  B  as  functions  of  x  .  The  first  two  terms  on  the  right 
hand  side  are  required  on  the  basis  of  the  results  for  ring  wings  at  incidence. 
The  necessity  of  the  additional  cosine  term  will  become  apparant  when  we  discuss 
the  propeller  tern  v^>  .  prom  the  Biot-savart  law  it  can  then  be  shown,  cf. 
Ref.  1  ,  that  (v  +v  *)  is  co-Dosed  of  an  axisymmetric  contribution  from  yQ  , 
plus  sin  0  and  cos  0  contributions  from  and  -37j/5f>  ,  or 

vv  “ ^  +osic  of 

where  x  a  x/r  ,  Axv  a  x-.xy  ,  X  a  c/2R  and 

3C  +  1/4  *  AXV[- Q^-Qjj+Qjj+Qjj]  /  4tt  (6) 

The  Q'-ts  aro  Legendre  functions  of  the  second  kind  and  order  3/2,  -1/2,  1/2, 
-1/2  respectively,  the  primes  denoting  differentiation  with  respect  to  their 
argument. 


*v> 


dx„  +  0  cos  0 


i 


JTb(xv)  dxy 


(5) 


CD  -  1  +  (A5v)2/2  (7) 

The  q'  terms  in  can  be  traced  to  the  bound  vortices  of  7^  and  the  other 

three  terms  to  the  corresponding  trailing  vortices  of  «d7j/30  . 

Tusning  to  the  remaining  shroud  contributions  to  v  due  to  the  unsteady 


part  7  ,  ve  have 

o 

s 

(8) 

Vy  »  0 

(9) 

v».  «  0 

(10) 
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This  follows  immediately  from  a  generalization  of  Eq.  (5)  for  the  expected  form 
of  7  . 

Finally,  we  take  up  the  various  propoller  contributions.  It  is  not  diffi¬ 
cult  to  prove  that  the  bound  contributions  vanish,  regardless  of  the  form  of  r 
and  , 

vr  «  0  (11) 

-  0  (12) 

and  that  the  trailing  contribution  v^/  is  independent  of  0  ,  so 

VP'  ”  vr^(*)  (W) 

This  effective  propeller  induced  camber  is  the  same  as  for  the  case  of  zero 
incidence . 

For  the  last  two  contributions,  we  assume  varies  sinusoidally  as  fit 
since  only  the  first  "harmonic"  is  present  in  either  the  case  of  the  duct  at 
incidence  or  of  the  propeller  at  incidence.  The  amplitude  and  the  phase  angle 
are  not  known  a  priori,  but  must  be  determined  by  solving  the  unsteady  propeller 
problem  explicitly.  However,  they  will  be  regarded  as  known  quantities  for  the 
moment.  From  Biot-savart  integration  then  we  find, 

v^/  ■  t.T  &  sin  9  +  XJ  &  cos  9  (14) 

Where  &  and  ^  are  functions  of  x  involving  complicated  double  integrals 
over  the  amplitude  and  phase  angle  of  .  Both  and  owe  their  pres¬ 

ence  to  the  helical  nature  of  the  propeller  vortex  system,  cf.  Ref.  2  ,  and 
tend  to  zero  in  the  limit  as  the  advance  ratio  of  the  propeller  becomes  Infinite. 
They  will  be  dealt  with  later. 

This  leaves  only  the  shed  terra  v» ,  from  the  propeller.  By  means  of  the 
result  for  v^ 

v^.  -  0  (15) 

That  is,  the  shed  vortex  system  may  be  interpreted  as  a  sequence  of  propellers 
of  incremental  loading,  displaced  both  axially  and  circumferentially  downstream 
of  the  propeller  plane,  each  one  of  which  is  governed  by  Eg.  (12)  . 
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1.4  g-Dccoroposltlon  of  Tangent-glow  Condition 


With  v  given  by  Eqs.  (2),  (3).  (5)  and  (8) -(15)  ,  the  tangent-flow  con¬ 
dition  of  Eg.  (1)  must  be  valid  for  like  terns  in  9  .  For  the  zeroth  terms, 
we  have 


+  v 


^0 


+  v  ,)/U 

r0 


(16) 


and  for  the  sin  6  and  cos  9  terms  respectively, 
■  X 
-X 


-  a  -  SF  =  J  (Axy)  A(xv)  dicv 


L 


<8  “  /  J^(Axv)  B(xv)  dxv 

-X 


(17) 


(18) 


where  is  the  kernel.  It  is  identical  to  our  -K^  of  Refs.  1  and  2  with 

mN  «•  1  and  u/flR_  •*  »  ,  which  has  been  related  to  Weissinger's  ring-wing 
7  P 

kernel' . 

We  see  now  that  Eg.  (16)  represents  the  same  ducted  propeller  configuration 
in  forward  flight  at  zero  incidence,  litis  equation  has  already  been  solved  for 
the  corresponding  shroud  vortex  and  source  distributions  in  Refs.  3  and  4  .  We 
therefore  turn  our  attention  to  the  "incidence"  equations  or  Eqs.  (17)  and  (18)  . 
As  we  have  said,  both  and  $  depend  on  the  amplitude  factor  and  phase 

angle  of  the  unsteady  part  of  the  propeller  circulation,  since  they  are  found 
to  depend  in  turn  on  the  unknowns  A  and  B  ,  these  equations  are  coupled 
through  $£  and  (S  . 
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CHAPTER  TWO 


SUPERPOSITION  PRINCIPLE  AS  A  FIRST  APPROXIMATION 

2.1  The  suporposltlon  Principle 

It  appears  physically  reasonable  that  as  a  first  approximation  A  and  B 
can  be  evaluated  by  using  a  simple  superposition  of  the  given  ducted  propeller 
at  zero  incidence  plus  a  cylindrical  duct  at  incidence  a  ,  oo?  Fig.  2  .  In 
terms  of  our  incidence  equations,  this  is  equivalent  to  neglecting  both  9 
and  entirely.  Thus  A  is  governed  by  the  "inclined-duct  equation", 

-  a  -  f  & (Ax;)  A(xv)  dxv  (19) 

and  B  is  identically  zero. 

Although  other  results  sure  available,  we  w^ll  outline  a  solution  to  Eq.  (19) 
based  upon  our  Legendre  function  formulation,  which  is  more  convenient  for  com¬ 
putation  and  for  subsequent  inclusion  of  the  propeller  interaction  terms  9 
and  (S  • 

2.2  Solution  for  the  Inclined  Duct 

We  introduce  the  Glauert  angular  variable  <p  ,  0  <  <t> «  it  ,  according  to  the 

relations, 

x  s  -  Xcos  <fi  ,  xv  e  -Xcos  0V  (20) 

and  seek  a  Glauert  series  solution  of  the  form 

00 

A  »  a  j  a^cot  -  +  ^  avsin  v0  1  (21) 

*•  v=l  -* 

which  satisfies  the  Kutta  condition  at  the  trailing  edge  and  has  the  proper 
square-root  singularity  at  the  leading  edge. 

As  in  the  case  of  the  ducted  propeller  at  zero  incidence,  the  kernel  be¬ 
haves  as  l/27rAxv  as  -*  0  ,  or 

CfC  a  1/27 TAXV  +  9t  (22) 

The  regular  part  91  is  equal  to  1/4  plus  an  antisymmetric  function  of  Axv 


ISITIOK  PRINCIPLE  FOR  QUOTED  PROPELLER  AT  ANGLE  OF  ATTACK 


which  approaches  1/4  as  Ax^  -*  ®  ,  see  Fig.  3  »  and  may  be  expanded  in  a 
double  Fourier  cosine  series  in  0,  0v  .  we  have  computed  the  first  eight 
Fourier  coefficients  in  each  series  numerically,  for  X  ■  0.25,  0.5,  0.75,  1.0  , 


FIGURE  3 

REGULAR  PART  OF  THE  KERNEL 


to  an  accuracy  of  ±1  in  the  fifth  significant  figure.  These  coefficients 
are  related  simply  to  Weiasinger's  b^  coefficients  and  agree  well  with  the 
corresponding  values  that  have  been  calculated  by  him  and  the  Bureau  Technique 
Zborowaki. 

With  thus  expressed,  Eqs.  (19)  -  (22)  produce  an  infinite  set  of  alge¬ 
braic  equations  in  the  unknown  Glauert  coefficients.  This  set  may  be  inverted 


*  The  authors  gratefully  acknowledge  numerical  checks  provided  through  private 
communication. 
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in  the  matrix  form 


(a)  -  -  [0]  (  1  *  0  ,  0  ,  ...  f  (23) 

where  (a)  is  a  column  vector  with  elements  aQ  ,  a^  ,  ....  and  the  inveroion 
matrix  [ 0  ]  is  an  infinite  square  matrix  with  elements  which  depend  on  X  . 
Retaining  only  aQ  through  ag  ,  we  have  inverted  the  truncated  system  numeri¬ 
cally  for  X  •  0.25,  0.5,  0.75*  1.0  and  have  tabulated  the  resulting  (a) 
vectors  in  Table  1  .  The  inversion  matrix  [ 0  ]  is  also  given,  in  Tables  2  -  5  , 
since  it  will  be  needed  later  in  the  analysis.  For  i  or  J  a  7  ,  the  elements 

0 .  .  may  be  approximated  by  their  asymptotic  values,  0 ,  .  —  2.0  and 
_*«j  , 

—  0.0  for  i  ^  J  . 

2.3  Shroud  Surface  Pressures 

The  steady  shroud  surface  pressure  coefficient  is  defined  in  the  conven¬ 
tional  form. 


(24) 


where  p  is  the  steady  static  pressure  on  the  shroud,  p  the  static  pressure 

00 

at  infinity,  and  p  the  fluid  mass  density.  It  follows  from  the  linearized 
Bernoulli  equation  that 


cp  «  -  2u/u  (25) 

Within  the  superposition  approximation,  this  may  be  rewritten  as  the  sum  of  the 
pressure  on  the  given  ducted  propeller  at  zero  incidence,  plus  the  pressure  on 
the  inclined  duct,  that  is. 


c  -  c  -  2u  (x,l-,0)  / u  (26) 

P  P0  7i 

with  *  DA  sin  0  .  Only  the  bound  vortices  on  the  duct  contribute  to  the 
inclined  duct  pressure  since  the  trailing  vortices  are  parallel  to  the  x-axis 
and  do  not  induce  any  axial  component  of  velocity. 

To  calculate  the  axial  velocity  induced  on  the  inclined  duct,  we  first 
consider  an  arbitrary  field  point  (x,r,0)  off  the  shroud  and  then  let  f  •*  1-  . 
In  the  limit,  this  reduces  to  the  local  "two-dimensional"  value  +H  UA  sin  8 
plus  a  "three-dimensional"  term.  The  three-dimensional  term  is  equal  to  U  sin  0 
times  the  Chordwise  integral  of  A(xv)  weighted  by  (Q^  +  -  2q£)/4it  , 

Eq.  (7)  giving  the  argument  of  the  Q's  . 
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It  la  convenient  to  express  this  result  in  matrix  form  for  use  in  conjunc¬ 
tion  with  {a}  from  Bq.  (23)  .  The  matrix  form  of  the  two-dimensional  part 
can  be  written  down  immediately.  For  the  matrix  form  of  the  three-dimensional 
part,  we  follow  Hough^.  Thus  Eq.  (26)  becomes 

CP  “  CP0  +  ®(±C*)T-  (a)  sin  6  (27) 

where  we  have  introduced  the  vectors, 

(<p}T  s  {  cot  -  ,  sin  0  ,  sin  20  ,  ...  )  (28) 

T 

(0)  r!  {  1  ,  -COS  0  ,  —COS  20  ,  ...  )  (29) 

Values  of  the  vectors  (q>)  for  x/c  =  -0.5,  -0.4,  ...  0.5  respectively  are 
given  by  the  rows  of  the  matrix  composed  of  COLUMNS  118,  120,  ...  130  of 
WORKSHEET  II  of  Ref.  4  }  similarly,  (0)  may  be  obtained  from  COLUMNS  205, 

207,  •••  217  .  We  have  calculated  the  elements  of  [8P]  numerically  to  within 
±1  in  the  fourth  significant  figure  and  tabulated  the  results  in  Tables  6-9 
for  X  »  0.25,  0.5,  0.75,  1.0  and  i,J  =  0  ,  ...  6  .  The  [&)(&)  vector  has 

also  been  calculated  and  is  included  in  Table  1  for  convenience. 

The  zero  incidence  contribution  to  the  pressure  coefficient  may  also  be 
written  in  similar  form,  consisting  of  the  two-dimensional  part  ±C(x)  Which 
gives  a  Jump  in  pressure  across  the  shroud  plus  a  three-dimensional  part  Which 
is  continuous  across  the  shroud  as  above,  see  Refs.  3  and  4  .  Like  a(x)  , 

c(x)  may  be  written  as  (<p)T(c)  Where  (c)  is  a  column  vector  of  the  olauert 

coefficients  for  the  vortex  distribution  on  the  given  ducted  propeller  at  zero 
incidence. 


2.4  Shroud  Forces  and  Moments 


With  the  pressure  coefficient  in  hand,  the  gross  forces  and  moments:  for 
the  shroud  can  be  determined.  Here  we  wi:  limit  ourselves  for  simplicity  to 
the  first-order  quantities,  namely  the  lift  and  its  associated  pitching  moment. 
By  symmetry  about  the  xz -plane,  there  is  no  side  force  or  yawing  moment. 


For  first-order  quantities,  we  need  only  the  net  loading 
positive  radially  outward.  Bq^  (27)  gives 


°P1 


cp]  •  “  2(<p)T{c)  -  2a(<p)T(a)  sin  6  (30) 

Which  has  the  same  chordwise  variation  as  c  ]  for  a  two-dimensional  airfoil. 
Therefore,  the  resulting  elemental  radial  force  and  moment  about  the  lending 


1? 


edge  on  the  strip  Rd 9  by  c  follow  immediately  from  thin  airfoil  theory 
Resolving  these  as  necessary  and  then  integrating  on  0  ,  we  find 


C/  s  - 5—0  *  -  7TX{2a0+  a,)a 

*  %pu  its 


sn  tx\ 

C„  b  - - — -  -  -x-  (2a0  + 2a1  -  ag)a 


m  JjpU2^3  2 


(31) 

(32) 


where  /  is  the  total  lift  on  the  shroud  and  tn  is  the  total  pitching 
moment,  positive  nose-up,  about  an  axis  parallel  to  the  y-axia  and  through 
the  leading  edge. 

As  we  would  expect,  c,  and  c.  are  the  aame  as  for  then  inclined  duct 

/  Tn 

alone.  We  hasten  to  point  out,  however,  that  Eq.  (32)  ia  correct  only  to  first 
order.  Practically,  there  are  second-order  terms  due  to  the  axial  forces  on 
the  shroud  Which  can  give  substantial  contributions  in  addition,  cf.  Kriebeli3. 


2.5  Propeller  Forces  and  Moments 

According  to  the  superposition  model  the  unsteady  loading  on  the  propeller 
is  zero.  Consequently,  the  only  forces  and  moments  we  have  are  the  usual  thrust 
and  torque  as  found  in  the  case  of  zero  incidence,  see  Ref.  4  . 
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0 
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E 
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(a) 


ss 

0.25 

0.50 

0.75 

0 

-1.4133 

-1.0977 

-0.9139 

-0.7955 

i 

0.0985 

0.2387 

0.3456 

0.4154 

2 

0.0040 

0.0180 

0.0409 

0.0686 

3 

-0.0003 

-0.0015 

-0.0021 

-0.0001 

4 

0.0000 

-0.0003 

-0.0012 

-0.0024 

5 

0.0000 

0.0000 

-0.0001 

-0.0005 

6 

0.0000 

0.0000 

0.0000 

0.0000 

[^JW 


IS 

0.25 

0.50 

0.75 

1.00 

0 

-0.2125 

-0.1746 

-0.1356 

-0.1099 

1 

0.1627 

0.2161 

0.2246 

0.2159 

2 

0.0089 

0.0349 

0.0636 

0.0866 

3 

0.0006 

0.0041 

0.0114 

0.0212 

4 

0.0000 

0.0000 

0.0006 

0.0024 

5 

0.0000 

0.0000 

0.0000 

0.0000 

6 

0.0000 

0.0000 

0.0000 

-0.0001 

TABLE  1 

{a}  AHD  [<?\ (a)  VECTORS  FOR  SUPERPOSITION  SOLUTION 
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f 

it 


•tm 

t 

> 

1.4133 

-0.2675 

-0.0256 

•  0.0004 

0.0001 

0.0000 

0.0000 

-0.0985 

1.9462 

0.0018 

0.0028 

C.0000 

0.C 

0.0000 

1 

A 

-0.0040 

0.0008 

1.9965 

0.0000 

0.0007 

0.0000 

0.0000 

t 

\ 

i 

0.0003 

0.0009 

0.0000 

1.9989 

0.0000 

'.0003 

0.0000 

*■"* 

0.0000 

0.0000 

0.0004 

0.0000 

1.9995 

0.0000 

0.0002 

i 

i . 

0.0000 

0.0000 

0.0000 

0.0002 

0.0000 

1.9997 

0.0000 

r 

1 

0.0000 

0.0000 

0.0000 

0.0000 

0.0001 

0.0000 

1.9998 

1 


k 


TABLE  2 

r 

INVERSION  MATRIX  [  (?)  ] 

FOR  X  a 

0.25 

r 

i. 

r 

1.0977 

-0.3811 

-0.0674 

-0.003? 

0.0004 

O.OOOl 

0.0000 

r 

-0.2387 

1.8508 

0.0147 

0.0156 

-0.0001 

-0.0003 

0.0000 

1 . 

-0.01&V, 

0.0062 

1.9812 

0.0001 

0.0036 

0.0000 

-0.0001 

r 

0.0015 

0.0044 

-0.0001 

1.9949 

0.0000 

0.00.  r 

0.0000 

v 

* 

i 

0.0003 

-0.0001 

0.0018 

0.0000 

1.9976 

0.0000 

0.0009 

[ 

0.0000 

-0.0001 

0.0000 

0.0009 

0.0000 

1.9986 

0.0000 

[ 

r 

0.0000 

o.ouoo 

0.0000 

0.0000 

0.0006 

0.0000 

1.9991 

[ 

r 

TABLE  3 

L 

INVERSION  MATRIX  [ 0  ] 

FOR  X  = 

0.50 

r 

L 
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r 

r 


r 

0.9139 

-0.4268 

-0.1072 

-0.0101 

0.0006 

0.0003 

0.0000 

r 

-0.3456 

1.7469 

0.0405 

0.0412 

-0.0002 

-0.0012 

0.0000 

l 

-0.0409 

0.0191 

1.9497 

0.0005  . 

0.0106 

0.0000 

-0.0004 

r 

1 

0.0021 

0.0115 

-0.0002 

1.9858 

0.0000 

0.0041 

0.0000 

r— , 

1 

0.0012 

-0.0006 

0.0052 

0.0000 

1.9939 

0.0000 

0.0022 

1 

0.0001 

-0.0003 

0.0000 

0.0025 

0.0000 

1.9966 

0.0000 

r 

0.0000 

0.0000 

-0.0001 

0.0000 

0.0003 

0.0000 

1.9990 

[ 

TABLE  4 

r 

INVERSION  MATRIX  [0] 

FOR  X  » 

0.75 

r 


r 

- 

1. 

0.7955 

-0.4418 

-0.1406 

-0.0206 

-0.0002 

0.0007 

0.0002 

r~’ 

1 

-0  .4154 

1.6448 

0.0734 

0.0767 

0.0001 

-0.0025 

-0.0001 

L  . 

r 

-0.0686 

0.0381 

1.9040 

0.0018 

0.0230 

-0.0001 

-0.0011 

[ 

0.0001 

0.0219 

0.0000 

1.9696 

0.0000 

0.0089 

0.0000 

r 

0.0024 

-0.0013 

0.0111 

-0.0001 

1.9872 

0.0000 

0.0044 

r 

0.0005 

-0.0007 

»C« '001 

0.0053 

0.0000 

1.9932 

0.0000 

[ 

0.0000 

0.0000 

X03 

0.0000 

0.0030 

0.0000 

1.9957 

c 


TABLE  5 

*  INVERSION  MATRIX  [0]  for  \  -  1.00 

c 

0 

- 
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0.1557 

0.0768 

0 

0.0011 

0 

0.0000 

0 

-0,1153 

0 

-0.0575 

0 

-0.0002 

0 

0.0000 

-0.0043 

0.0287 

0 

-0.0309 

0 

0.0000 

0 

-0.0004 

0 

0.0206 

0 

-0.0208 

0 

0.0000 

0.0000 

0.0000 

Q 

0.0157 

0 

-0.0158 

0 

0.0000 

0 

0.0000 

0 

0.0125 

0 

-0.0125 

0.0000 

0.0000 

0 

-0.0002 

0 

0.0106 

0 

TABLE  6 

MATRIX  FOR 

X  “  0.25 

0.1772 

0.0833 

0 

0.0053 

0 

0.0000 

0 

-0.1985 

0 

-0.0977 

0 

-0.0016 

0 

0.0000 

-0.0211 

0.0489 

0 

-0.0590 

0 

-0.0C04 

0 

-0.0031 

0 

0.0393 

0 

-0.0408 

0 

-0.0001 

0.0000 

0.0004 

0 

0.0306 

0 

-0.0310 

0 

0.0000 

0 

0.0001 

0 

0.0248 

0 

-C.0249 

0.0000 

0.0000 

0 

0.0000 

0 

0.0208 

0 

TABLE  7 

[^]  MATRIX  FOR  X  »  0.50 


0.1775 

0.0771 

0 

0.0116 

0 

0.0001 

0 

-0.2512 

0 

-0.1206 

0 

-0.0050 

0 

0.0000 

-0.0466 

0.0603 

0 

-0.0822 

0 

-0.0014 

0 

-0.0099 

0 

0.0548 

0 

-0.0593 

0 

-0.0005 

-0.0003 

0.0012 

0 

0.0444 

0 

-0.0456 

0 

0.0000 

0 

0.0005 

0 

0.0365 

0 

-0.0370 

0.0000 

0.0000 

0 

0.0002 

0 

0.0308 

0 

TABLE  S 

\.8>\  MATRIX  FOR  / 

»  0.75 

0.1737 

0.0682 

0 

0.0183 

0 

0.0004 

0 

-0.2827 

0 

-0.1310 

0 

-0.0104 

0 

0.0000 

-0.0747 

0.0655 

0 

-0.0995 

0 

-0.0034 

0 

-0.0207 

0 

0.0663 

0 

-0.0754 

0 

-0.0013 

-0.0016 

0.0026 

0 

0.0565 

0 

-0.0593 

0 

0.000 r 

0 

0.0013 

0 

0.0475 

0 

-0.0488 

0.0001 

0.0000 

0 

0.0006 

0 

0.0404 

0 

TABLE  9 

[&J 

MATRIX  FOR 

A  *  1.00 
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CHAPTER  THREE 


INCLUSION  OF  SHROUD-PROPELLER  COUPLING 

3.1  The  "Interaction  Model" 

Although  the  simple  superposition  hypothesis  appears  to  be  reasonable  on 
physical  grounds,  it  is  not  easily  justified  by  the  usual  order  of  magnitude 
type  reasoning.  Therefore  we  will  investigate  the  effect  of  the  shroud-propeller 
coupling,  manifested  by  the  presence  of  S  and  *§  in  the  full  integral  equa¬ 
tions,  Eqs .  (17)  and  (18)  . 

The  only  way  this  can  be  carried  out  practically  is  by  an  iteration  pro¬ 
cedure.  That  is,  with  A  and  B  (=0)  given  by  the  superposition  model  as  a 
first  approximation,  we  calculate  the  corresponding  6-variation  in  inflow,  or 
periodic  gust,  seen  by  each  blade  olement  as  the  propeller  rotates,  we  then 
solve  the  associated  unsteady  propeller  problem  for  V ^  .  With  ,  S  and 
rS  follow  and  Eqs.  (17)  and  (18)  become  equivalent  to  the  equations  of  an 
inclined  duct  with  non-axiaymmetric  camber,  consequently,  these  equations  can 
be  inverted  in  the  same  fashion  and  so,  a  second  approximation  found  to  A  and 
B  .  The  procedure  may  in  turn  be  repeated  successively. 

Since  S  and  turn  out  to  be  fairly  small  compared  to  a  ,  we  limit 

the  calculations  for  the  present  purposes  to  the  second  approximation.  In 
addition,  we  assume  f .  is  uniform  over  the  blade  and  determine  the  amplitude 

A  " 

and  phase  of  by  applying  classical  unsteady  two-dimensional  strip  theory 
at  a  representative  radius. 

3.2  Determination  of  unsteady  inflow  and  Propeller  circulation 

We  fix  ourselves  for  the  moment  to  the  fth  blade  at  the  nominal  point 
(xp,r,0£)  ,  Whore 

6j  h  -  flt  +  2ir£Al  (33) 

so  that  the  0th  blade  coincides  with  the  direction  of  the  y-axls  at  t  =>  0  . 
With  cos  a  <3  1  ,  then  the  blade  element  centered  about  this  point  sees  an 
effective  free  stream 


Ql  b  V  u2  +  n2r2  (3^) 

along  the  local  helix  angle,  together  with  unsteady  axial  and  tangential  com¬ 
ponents  of  inflow. 

Within  our  approximations,  the  unsteady  axial  component  is  due  to  the  shroud 
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bound  vortices  and  the  unsteady  tangential  component,  to  the  incidence  crossflow 
and  the  shroud  trailing  vortices.  We  resolve  these  components  parallel  and 
perpendicular  to  Ql  .  The  effect  of  the  parallel  or  surging  component  is  small 
and  may  be  neglected.  For  the  perpendicular  or  unsteady  downwaah  component, 
we  find 


to  «  u 


Jlr 


yl<w 


-|ua 


cos  0.  +  w.  , 

t  y1 


)w 


(35) 


Where  to  is  taken  positive  if  directed  "downward".  Analogous  to  Eg.  (5)  ,  It 
can  be  shown  that  the  axial  velocity  induced  by  the  bound  vortices  of  the  duct 
loading  7^^  n  UA  sin  0  is  proportional  to  U  a  sin  0  ,  and  that  the  tangential 
velocity  Induced  by  the  associated  trailing  vortices  is  proportional  to 
Uacos  0  .  Setting  0  °  0^  and  combining  the  latter  term  with  the  crossflow 
term,  we  can  rewrite  Eg.  (35)  as 


(36) 


Where  J  =  u/HR  is  the  advance  ratio  of  the  propeller,  u  s  r /r  and 
P  P 


(37) 


dx_. 


(38) 


with  A(xv)/a  given  by  Egs.  (21)  and  Table  1  .  The  argument  3  and  modulus 
jc  of  the  Heuman  lambda  function  A0(3»x)  axe 


3  ■=  sin 


-1 


(A*v)‘ 


<AXV)2  +  (l-?)2 


4r 

V  (AiL)2  +  (1+r) 


(39) 


and  the  argument  of  the  Legendre  functions  is 


ffl  “  1  +  [(Axy)2  +  (l-r)2]/2r 


(*0) 


For  r  =  0.7  ,  the  constants  1^  ,  1^  have  been  evaluated  by  numerical  inte¬ 
gration  for  several  values  of  X  and  xp/c  to  an  accuracy  of  +0.0001  ,  and 
are  tabulated  in  Table  10  . 
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Xp/c 

X  \ 

-0.25 

-0.125 

0.0 

0.125 

0.25 

0.25 

0.5506 

0.5201 

0.4696 

0.4088 

0.3449 

0.50 

0.6182 

0.5113 

0.4021 

0.3046 

0.2207 

0.75 

0.5445 

0.3938 

0.2753 

0.1867 

0.1200 

1.00 

0.4488 

0.2863 

0.1789 

0.1090 

0.0624 

\a/c 

X 

-0.25 

-0.125 

0.0 

0.125 

0.25 

0.25 

0.8195 

0.8568 

0.8985 

0.9413 

0.9830 

0.50 

0.8418 

0.9298 

1.0119 

1.0846 

1.1474 

0.75 

0.9128 

1.0327 

1.1308 

1.2090 

1.2708 

1.00 

0.9894 

1.1248 

1.2248 

1.2978 

1.3511 

TABLE  XO 


nr-OCED  VELOCITY  FACTORS 


\  ' 


1^  FOR  r  =  0.7 


The  chordwise  point  we  chose  on  the  blade  element  was  arbitrary.  If  we 
had  taken  another  point  at  a  distance  g  further  along  the  chord  in  the  tfl- 
direction,  to  would  still  be  given  by  Eq.  (36)  but  with  9  ^  replaced  by 
9 1  plus  the  increment  .  Recalling  Eq.  (33)  then,  we  see  this  resulting 

downwaah  is  equivalent  to  the  sinusoidal  gust  problem  solved  by  N.  Kemp1-’  and 
others.  Taking  d  as  the  chord  length  of  the  element  and  adapting  his  solution, 
we  get 


»  oOd  ocoa(9^-0) 


(41) 
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with  the  amplitude  factor  a  and  phase  angle  0  given  by 


-  [vr-  £  vf  +  [v1  + |  vr]‘ 


i>  -  tan 


-1 


JuV*1  +  ?Vrl 
|_JnVr  -  f^j 


j*  4 

The  factors  a  and  a  depend  only  on  the  reduced  frequency 


k  b  ( ld/2 <%£ 


and  in  terms  of  the  Bessel  and  Hankel  functions  are  defined  by 


.  ..  J,  +  iJn 

a  b  ar  +  ia1  s  (2e"1V7rk)  'J 


+  IHq2^ 


These  are  tabulated  in  Table  11  . 


(42) 

(43) 


(44) 


(45) 


3*3  Calculation  of  Shroud-propeller  Coupling 

For  uniform  f j  ,  we  are  concerned  only  with  concentrated  helical  vortices 
*  trailing  from  each  of  the  N  blade  tips.  If  r  represents  the  parametric 
time  for  a  vortex  element  to  convect  downstream,  its  strength  at  any  point  is 
atJdo  cos(6^  -  ■$  +  ftr)  and  varies  periodically  along  the  length  of  the  helix. 

The  tota..  radial,  wash,  steady  and  unsteady,  induced  on  the  shroud  by  these 
vortices  may  then  bo  expressed  according  to  the  Biot-Savart  law  in  the  same 
fashion  as  for  a  helical  vortex  of  constant  strength  except  for  the  above 
cosine  dependence  Which  now  appears  in  the  t-J ntegration.  To  obtain  the  steady 
part,  wa  integrate  in  turn  with  respect  to  over  the  period  2-n  and  divide 

by  2ir  .  If  the  order  of  integration  is  inverted,  the  flt -integration  may  be 
performed  explicitly  in  terms  of  Legendre  functions.  The  result  is  independent 
of  the  blade  index  4  ,  so  that  the  summation  over  the  blade  number  merely 
produces  a  factor  N  .  It  may  therefore  be  shown  that 


vp  “  -  2HJL-Z?  | cosde  f  A^pT(Q4+  QJj)d^  +  JjisinAtf  f  (Q^-Qjj)^  j  (46) 

8w2Rp  (  Jq  Jq  ) 

with  AXpT  s  (x - Xp  - ^.t)  ,  ?sj)T,  fifl  3  (9-0)  and  the  argument  of  the  Q*s 
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k 

0.0000 

0.0400 

0.0800 

0.1200 

0.1600 

0.2000 

0.2400 

0.2800 

0.3200 

O.36OO 

0.4000 

0.4400 

0.4800 

0.5200 

0.6000 

0.7000 

0.8000 

0.9000 

1.0000 


1.0000 

0.9153 

0.8273 

0.7484 

0.6796 

0.6199 

0.5677 

0.5218 

0.4811 

0.4447 

0.4119 

0.3821 

0.3548 

0.3297 

0.2849 

0.2368 

0.1954 

0.1591 

0.1269 


0.0000 

-0.1697 

-0.2553 

-0.3065 

-0.3382 

-0.3581 

-0.3705 

-0.3779 

-0.3820 

-0.3836 

-0.3836 

-0.3823 

-0.3801 

-0.3772 

-0.3699 

-0.3590 

-0.3467 

-0.3335 

-0.3196 


0.0000 

-0.1145 

-0.1546 

-O.I678 

-0.1675 

-0.1596 

-0.1476 

-0.1333 

-0.1177 

-0.1014 

-0.0849 

-0.0684 

-0.0521 

-0.0360 

-0.0049 

0.0318 

0.0659 

0.0972 

0.1259 


br 

1.0000 

0.9240 

0.8526 

0.7926 

0.7429 

0.7016 

0.6667 

O.6369 

0.6110 

0.5882 

0.5679 

0.5495 

0.5326 

0.5169 

0.4884 

0.4561 

0.4260 

0.3971 

0.3686 


TABLE  11 

UNSTEADY  AIRFOIL  THEORY  COEFFICIENTS  a(k)  AND  b(k) 
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is  given  by 


a  s  X  +  [(AxpT)2  +  (i-n)2]/2p.  (47) 

The  first  integral  is  an  exact  differential  and  can  be  evaluated  immediately. 
Integration  of  the  second  integral,  however,  is  quite  difficult.  After  con¬ 
siderable  manipulation,  we  can  find  an  approximate  result.  The  approximation 
is  in  the  sense  that  p  is  close  to  unity,  which  is  true  in  practice.  Thus, 
we  have 

I 

Sf  «  -  —  |  ^(Qjj+  QJj)  sin'>  +  AiCp(Q^+Q_^)]  COSi>  |  (48) 

where  co  is  now  given  by  Eq.  (47)  vith  r  ”  0  .  The  function  is  identical 
to  S  ,  with  sin  0  and  cos  ■a  replaced  by  cos  and  -  sin  <7  respectively. 

Before  evaluating  S'  and  8  one  point  should  be  made  regarding  their 
form.  Notice  that  both  terms  are  proportional  to  H  .  It  would  therefore  seem 
that  for  a  given  value  of  the  steady  part  of  the  propeller  thrust,  they  could 
be  made  arbitrarily  large  —  and  hence  also  tht*  corresponding  shroud  loading  —  by 
increasing  the  blade  number.  Actually,  as  N  becomes  large,  the  effect  of 
blade  interference  would  come  into  play  and  decrease  the  amplitude  of  the  un¬ 
steady  propeller  circulation. 

To  illustrate  the  quantitative  and  qualitative  importance  of  the  shroud- 
propeller  coupling,  consider  the  typical  set  of  parameters  N  **  3  ,  U  “  0.95  » 


FIGURE  4 


d/Rp  ■  0.2  ,  Xp/c  »  -0.25  for  three  cases*  (i)  X  ■  0.5  *  J  *  0.5  »  (ii)  X  ® 
0.5  ,  J  «  0.25  7  (iii)  X  ■  1.0  ,  J  ■  0.5  .  we  have  calculated  S'  and  8 
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and  plotted  them  against  the  Giauevt  variable  In  Figs.  4  and  5  . 


«/ a 


FIGURE  5 


Since  their  Importance  is  determined  by  their  magnitude  relative  to  unity, 
we  see  that  the  shroud-propeller  coupling  is  in  fact  large  enough  to  account 
for  a  significant  improvement  over  the  simple  superposition  model,  but  not 
large  enough  to  discredit  the  superposition  model  as  a  reasonable  first  approxi¬ 
mation. 

Parametrically,  comparison  of  these  cases  indicates  an  Increase  in  the 
magnitude  of  the  coupling  terms  as  J  decreases  and/or  \  increases. 

3.4  The  Interaction  Solution 

The  second  approximation  or  interaction  solution,  i.e.  including  the  effect 
of  the  shroud-propeller  coupling,  may  be  constructed  easily  using  the  inclined- 

C 

duct  results  or  weissinger’s  general  formulation  .  We  take 

00 

A  »  a  1"  «gCot  |  ^  ayain  1  (49) 

v®l 

[co 

bQcot  |  ^  bvsin  1  (50) 

v-1  J 

and  expand  the  coupling  terms  as 

00 

&/<*  3  /o  "  /vcos  v*  (51) 

v-1 
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(52) 
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5  9q  ~  £  PvCOB  v<t> 
v**l 

Then,  it  in  not  difficult  to  show  that,  analogous  to  Eq.  (23)  , 

(a)  -  -  [0]  {  (1+  /Q)  ,  /j_  ,  /2  ...  }T  (53) 

P>)  a  -  (^]  {  q  ,  ^ •  ^2  '  ***  J  (54) 

Where  the  /v*3  and  yv'e  must  be  computed  by  numerical  Fourier  analysis. 
3.5  shroud  surface  Pressures 

Direct  extension  of  the  matrix  form  for  the  shroud  surface  pressures  from 
Eq.  (27)  yields, 

1  cv  "  co  +  «(+{(p]T-  {0)T[^J)  ((a)sin  9+  (b)cos  e)  (55) 

y  “ 0 

where  we  recall  that  the  values  of  [SP\  have  already  been  presented  in 
Tables  6-9  * 


3.6  shroud  Forces  and  Moments 

Consistent  with  the  superposition  model,  we  will  consider  only  first-order 
contributions  to  the  forces  and  moments. 

If  we  write  down  the  not  loading  cp]  ,  we  have  three  terms,  i.e.  an 
axisymmetric  term,  a  sin  8  term  and  a  cos  8  term.  The  axi symmetric  term  is 
identically  the  same  as  for  the  superposition  nodel  and  again  does  not  contribute 
to  any  gross  quantity.  The  sin  8  term  is  also  the  same  in  form  and  gives  a 
lift  and  associated  pitching  moment  as  expressed  by  Kqs.  (31)  and  (32)  ,  only 
the  av*s  of  Eq.  (23)  are  replaced  by  the  av's  of  Eq.  (53)  .  The  cos  8 
term  is  completely  new  and  has  not  appeared  in  other  analyses.  It  is  symmetric 
about  the  xy-plane  and  causes  a  side  force  A  and  yawing  moment  ta  which  by 
analogy  with  the  lift  force  and  pitching  moment  juts. 


c  a  - 

‘  •  V 


‘jpu2 


*  -  nX(2aQ+  a^Ja 


m  nX 

c  B  -  ■ - 

n  o  3  2 


Jjpu2 


7«- 


(2b0+2b1-b^)a 


(56) 

(57) 
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The  side  force  is  positive  in  the  positive  y -direction  and  the  yawing  moment, 
taken  about  an  axis  parallel  to  the  z-axie.  and  through  the  leading  edge,  is 
positive  in  the  nose-right  direction. 

3-7  Propeller  Forces  and  Moments 

For  the  superposition  model,  we  had  only  the  usual  thrust  and  torque  on  the 
propeller.  For  the  interaction  model,  though,  we  must  expect  additional  forces 
and  moments  corresponding  to  . 

To  compute  these  quantities,  the  local  unsteady  "lift"  is  needed  for  each 
blade  element.  From  Kemp1-’  again,  this  is  given  by  the  so-called  Kutta- 
Joukowski  lift,  times  the  factor  eiJ</[J0(k)-  iJ^(k) '  .  If  the  latter  is  com¬ 
bined  with  the  circulation,  however,  the  resulting  unsteady  lift  per  unit  length 
may  be  reduced  to  the  same  form  cfiUv  .  with  ar,i  in  a  and  0  of  Eqs.  (42) 

J*  1  W 

and  (43 )  replaced  by  b  '  where 

eik 

b  s  1“  +  ib1  s  a  •  -  (58) 

Jo-iJi 

The  quantity  b(k)  is  identical  to  "S(m)"  tabulated  by  Kemp1-*  and  partially 
reproduced  here  for  convenience  back  in  Table  11  .  Of  course,  b  ■*  a  as 
k  -  0  . 

If  we  resolve  the  above  incremental  lift  into  its  y-  and  z -components, 
integrate  ever  the  blade  radius,  and  sum  over  the  N  blades,  we  find  that  the 
associated  total  propeller  lift  and  side  force  are 


N-l 

L  +  L  -  apl^Rpdoj,  cos (0^  —  ob)  cos  0 ^  (59) 

1=0 

N-l 

S  +  S  =*  -  apU^RpdOk  coc(0^-db)  sin  0^  (60) 

f,°0 

m  j  4 

where  0^  s  a(b  '  )  and  5  -i(b')  .  Expressing  the  sines  and  cosines  in 

exponential  form,  the  steady  portions  of  these  summations  may  be  evaluated  ex¬ 
plicitly  as  Jj  N  cos  djj  and  >j  N  sin  i>b  respectively.  It  follows  that  the 
steady  lift  L  and  side  force  S  are 


[ 

[ 


CL  £ 


ijpU2 


7TR 


ap,2N(Jbd 

7TR_ 


cos 


(61) 
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(62) 


CS  3 


HpU2™2 


o 

an  Ntjj^d 

rrR 


sin  i>v 


The  resultant  force  vector  is  oriented  at  the  angle  0  =  (i>b  +  m/2) 
we  can  obtain  the  p itching  moment  M  and  yawing  moment  N  as 


M 


M  tyl^irR3 


a^Ni^d 

— — — —  sin  i). 
37TJR  b 


similarly, 


(63) 


C  S 


N  “I*  Nobd 


N  Jjpt^-rrR3  37rJRp 


cos  l>. 


(64) 


with  the  nose-up  and  nose-right  directions  taken  positive  as  for  the  shroud. 
Again  the  resultant  moment  vector  is  oriented  at  the  angle  6  *»  ( i>b  +  7t/2  )  . 

These  additional  forces  and  moments  are  closely  related  to  those  obtained 
by  programmed  blade  pitch  variation  as  studied  by  W.P.A.  Joosen  et  al1^  in 
connection  with  the  control  of  tandem-propeller  submarines. 
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NUMERICAL  RESULTS  AN2  COMPARISON  WITH  EXPERIMENT 


4.1  Cages  Chosen 

For  comparison  with  experiment,  we  have  chosen  Cusss  "2",  "3"  and  "4"  of 
Ref.  17  where 


Case 

\ 

(A 

J 

i»2 11 

0.345 

0.970 

0.355 

”3" 

0.500 

0.956 

0.344 

••Jim 

0.724 

0.936 

0.406 

with  N  *»  3  ,  xp/c  "  -0-219  and  d/Rp  °  0.213  in  each  case.  The  values  of 
J  and  d/Rp  ,  not  given  therein,  were  obtained  by  private  communication  with 
the  author. 

4.2  Theoretical  Results 

Let  us  outline  the  calculations  for  one  case,  say  case  "3"  .  For  tho 
superposition  model,  the  (a)  vector  is  given  directly  up  to  v  *■  6  by  Table  1  . 
The  surface  pressure  coefficient  is  found  then  from  Eg.  (27)  ,  and  the  lift  and 
moment  from  Eqs.  (31)  and  (32)  . 

For  the  solution  according  to  the  interaction  model,  we  first  compute  the 
reduced  frequency  k  -  0.132  from  Eq.  (44)  and  hence  obtain  ar  ®  0.727  and 
a*  »  -0-318  from  Table  11  .  From  Table  10  ,  the  inflow  factors  are  »  0.592 
and  ■  0.864  .  Eqs.  (42)  and  (43)  then  give  the  amplitude  factor  a  *  1.678 
and  phase  angle  <t  **  31-9°  -  With  these,  S’/a  and  *8/0,  may  be  calculated 
over  the  shroud  chord  from  Eq.  (48)  ,  and  f  ^  =  -0.1277  ,  0.0095,  -0.0224, 
-0.0184,  -0.0020,  0.0081,  0.0066  and  ■  -0.1123,  0.0575,  -0.0277,  -0.0313, 
-0.0059,  0.0120,  0.0111  obtained  by  numerical  Fourier  analysis,  see  Eqs.  (51) 
and  (52)  .  The  Glauert  coefficient  vectors  (a)  and  (b)  follow  from  Eqs.  (53) 
and  (54)  ,  and  hence  the  angle  of  attack  portion  of  the  net  surface  pressure 
coefficient  from  Eq.  (55)  •  With  (a)  and  (b)  ,  we  may  then  evaluate  the 
gross  duct  quantities  easily  from  Eqs.  (31)  ,  (32)  ,  (56)  and  (57)  «  For  the 
propeller,  however,  we  must  first  determine  ■  1.624  and  «*  43.4®  from 
Eqs.  (42)  and  (43)  ,  with  ar  and  a*  replaced  by  br  «  0.778  and  b*  ■  -0.168 
from  Table  11  .  In  turn,  the  steady  forces  and  moments  on  the  propeller,  which 
arise  due  to  the  unsteady  portion  of  the  blade  loading,  are  found  from 
Eqs.  (61)-  (64)  . 
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All  of  tbeso  results  for  tbs  Interaction  model  are  compared  with  the  re¬ 
sults  of  tha  superposition  mod-  .1  in  the  following  tables  and  Fig.  6  .  The 
overall  effect  oi  the  shroud -propeller  interaction  may  be  seen  clearly.  From 


Superposition  Interaction 

Model  Model 


V 

(a] 

0>) 

(a) 

(b) 

0 

-1.0977 

0 

-0.9555 

0.1432 

1 

0.2387 

0 

0.1912 

-0.1323 

2 

0.0180 

0 

0.0601 

0.0526 

3 

-0.0015 

0 

0.0354 

0.0623 

4 

0.0003 

0 

0.0038 

0.0119 

5 

0.0000 

0 

-0.0161 

-0.0240 

6 

0.0000 

0 

-0.0132 

-0.0222 

the  two  viluos  of  dc ^  /da  ,  we  find  that  it  causes  a  reduction  of  about  10#  in 
the  shroud  lift  or,  equivalently,  in  the  "effective”  incidence.  This,  however, 
is  approximated ' *  compensated  for  by  the  introduction  of  a  propeller  lift,  see 
dc  L  /da  ,  so  that  the  combined  lift  of  the  shroud  and  propeller  remains  about 
the  same.  This  compensatory  effect  ie  not  present,  though,  in  tho  case  of  the 
side  force  generated  by  the  interaction.  There  tho  shroud  and  propeller  side 
forces  are  additive  and  produce  a  resultant  force  equal  in  magnitude  to  about 
15#  of  the  lift,  see  dC^ /da  and  dCg/da  .  Regarding  the  moments ,  we  see  that 
the  interaction  increases  the  pitching  moment  by  about  20#  and  also  introduces 
a  yawing  moment  about  30#  as  big. 


Quantity 

Superposition 

Model 

Interaction 

Model 

dc  ^  /da 

3.074 

2.701 

dc  .  /da 

a 

0 

-0.242 

drVda 

WWW 

-0.682 

-0.624 

dc„  /da 

ww 

0 

0.013 

dc  L  /da 

0 

0.219 

dc  s  /da 

0 

-0.207 

dcM  /da 

0 

-0.192 

dcN/da 

0 

0.203 

30 
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FIGURE  6 

ANGLE  OF  ATTACK  CONTRIBUTION  TO  NET  SHROUD  PRESSURE  COEFFICIENT 
COMPARISON  OF  MODELS  FOR  CASE  "3” 
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We  have  omitted  the  corresponding  calculations  for  Cases  M2"  and  "4M  for 
brevity.  They  are  entirely  similar  to  Case  "3".  both  quantitatively  and 
qualitatively. 


4.3  Comparison  with  Experiment 

The  total  lift  3?  on  the  ducted  propeller  in  general  is  equal  to  the  sum 
of  the  duct  lift,  the  propeller  lift  and  the  component  of  the  total  thrust 
normal  to  the  stream.  Putting  this  in  coefficient  form  and  taking  d(  )/da 
of  the  result  gives 


d  d  d 

da  CSe  m  V  "  da  C/  +  da  CL 


(65) 


Experimental  values  for  the  terms  on  the  left  hand  side  (LHs)  are  given  in 
Kef.  17  .  For  comparison,  the  terms  on  the  right  hand  side  (RHS)  have  been 
computed  as  outlined.  These  results,  based  upon  the  interaction  model,  are 
presented  in  the  following  table. 


Case 

as  cse 

c ff" 

a  „ 
da  c/ 

9_  c 
da  UL 

LHS 

RHS 

Mg’* 

3.04 

0.2 

2.30 

0.23 

2.84 

2.53 

"3" 

3-34 

0.2 

2.70 

0.22 

3-14 

2.92 

3.61 

0.2 

3.11 

0.23 

3.41 

3.34 

The  corresponding  values  of  the  RHS  for  the  superposition  model  are  2.57*  3*07 
and  3*49#  Which  differ  very  little  from  the  values  predicted  by  the  interaction 
model.  Where  the  two  models  do  differ  significantly  though,  is  in  the  predic¬ 
tion  of  the  side  force  and  the  breakdown  of  the  total  lift  between  the  duct  and 
propeller.  Unfortunately,  these  new  forces  and  moments  are  not  reported  in 
Ref.  17  ,  nor  elsewhere  as  far  as  we  know,  and  must  await  further  tests  for 
substantiation . 
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A  theory  has  been  developed  for  the  steady  portion  of  the  aerodynamic  load¬ 
ing  on  a  ducted  propeller  In  forward  flight  at  angle  of  attack.  In  particular: 

Ka  have  shown  that  as  a  first  approximation  the  ducted  propeller  at 
angle  of  attack  may,  without  a  priori  assumption,  be  regarded  as  the 
superposition  of  the  given  ducted  propeller  at  zero  incidence  plus  a 
cylindrical  duct  of  the  same  aspect  ratio  but  zero  thickness  and 
camber  at  the  given  incidence.  Formulas  and  tables  are  presented  for 
the  convenient  calculation  of  the  cylindrical -duct  portion,  the  ducted 
propeller  at  zero  incidence  having  been  treated  in  previous  reports. 

This  approximation  has  also  been  refined  to  include  the  effect  of 
the  cyclic  variation  of  blade  loading,  due  to  the  asymmetric  flow 
through  the  propeller  plane,  and  the  associated  interaction  with  the 
shroud.  The  results  are  reduced  to  a  tractable  form  and  provide  an 
improvement  over  the  superposition  model  that  is  significant,  though 
not  great  enough  to  invalidate  the  superposition  theory  as  a  first 
approximation . 

More  important,  the  improved  "interaction"  model  predicts  a  net  side 
force  not  accounted  for  at  all  by  the  simple  superposition  model. 

It  may  be  as  great  as  20$  or  more  of  the  lift;  approximately  half  on 
the  propeller  and  half  on  the  shroud.  A  decrease  in  shroud  lift  of 
around  10$  is  also  predicted  by  the  interaction  model,  although  it 
is  approximately  balanced  by  a  lift  force  on  the  propeller. 

Comparison  with  tests  shows  good  agreement  for  the  values  of  total 
ducted  propeller  lift.  Unfortunately,  neither  side  forcer  nor  pro¬ 
peller  lift  have  been  measured  experimentally  so  that  d4  .  verifi¬ 
cation  of  the  additional  loads  predicted  by  the  inter a.  -n  model  is 

not  possible  at  this  time. 

In  the  light  of  these  new  results,  we  recommend  that  certain  aspects  of  the 
overall  aerodynamic  design  of  VTOL  ducted-propeller  aircraft  be  re-examined  as 
well  as  the  detailed  mechanical  design  of  the  shroud  and  propeller  components. 
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